Modeling exchange bias microscopically 

U. Nowak, A. Misra and K. D. Usadel 
Theoretische Physik, Gerhard-Mercator-Universitat Duisburg, 
47048 Duisburg, Germany 
e-mail: uli@thp.uni-duisburg.de 

February 1, 2008 



Abstract 

Exchange bias is a horizontal shift of the hysteresis 
loop observed for a ferromagnetic layer in contact 
with an antiferromagnetic layer. Since exchange 
bias is related to the spin structure of the an- 
tiferromagnet, for its fundamental understanding 
a detailed knowledge of the physics of the anti- 
ferromagnetic layer is inevitable. A model is in- 
vestigated where domains are formed in the vol- 
ume of the AFM stabilized by dilution. These 
domains become frozen during the initial cooling 
procedure carrying a remanent net magnetization 
which causes and controls exchange bias. Varying 
the anisotropy of the antiferromagnet we find a 
nontrivial dependence of the exchange bias on the 
anisotropy of the antiferromagnet. 
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1 Introduction 

When a ferromagnet (FM) is in contact with an 
antiferromagnet (AFM) a shift of the hysteresis 
loop along the magnetic field axis can occur which 
is called exchange bias (EB). Usually, this shift is 
observed after cooling the entire system in an ex- 
ternal magnetic field below the Neel temperature 
T N of the AFM. Although this effect is well known 
since many years Q and is already intensively ex- 
ploited in magnetic devices its microscopic origin 
is still under debate. For a review of the experi- 
mental work see the recent article by Nogues and 
Schuller @. 



In the approach of Malozemoff || exchange bias 
is attributed to the formation of domain walls in 
the AFM perpendicular to the FM/AFM interface 
due to interface roughness. However, the forma- 
tion of domains in the AFM only due to interface 
roughness is unlikely to occur because the creation 
of the domain walls is energetically unfavorable. 

Koon considered a spin-flop coupling between 
FM and the compensated AFM as responsible for 
EB |j, but recently, Schulthess and Butler §,§ 
showed that spin-flop coupling alone cannot ac- 
count for this effect. Instead, in their model EB 
is only obtained if uncompensated AFM spins are 
assumed at the interface — their occurrence is not 
explained microscopically. 

In a previous Letter Miltenyi et al. showed 
that it is possible to strongly influence EB 
in Co/CoO bilayers by diluting the antiferro- 
magnetic CoO layer, i. e. by inserting non- 
magnetic substitutions (Coi-^Mg^O) or defects 
(Coi-yO) not at the FM/AFM interface, but 
rather throughout the volume part of the AFM. 
In these systems the observed EB is primarily not 
due to disorder or defects at the interface. Rather, 
the full antiferromagnetic layer must be involved 
and it was argued that in these systems EB has its 
origin in a domain state in the volume part of the 
AFM which triggers the spin arrangement and the 
FM/AFM exchange interaction at the interface. 
This domain state develops due to the dilution 
of the AFM: the domain walls pass preferentially 
through non-magnetic sites thus reducing consid- 
erably the energy necessary to create a wall. The 
domain state strongly depends on the dilution of 
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the AFM resulting in a strong dependence of EB 
on dilution. Since dilution favors the formation of 
domains it leads to an increase of the magnetiza- 
tion in the AF and thus to a strong increase of the 
EB upon dilution (see also |8| , where it was shown 
that it is possible to influence (to increase or even 
to revers) EB by a subsequent ion irradiation of 
the sample). 

In the same letter this picture was further sup- 
ported by Monte Carlo simulations. Later it was 
shown |9|[l0| that a variety of experimental facts 
associated with exchange bias can be explained 
within our model, like positive exchange bias after 
cooling in strong magnetic fields, the temperature 
dependence of exchange bias, especially the rela- 
tion between the so-called blocking temperature 
and the Neel temperature, and the training effect, 
among others. In these studies the AFM CoO 
investigated experimentally was due to its rather 
strong uniaxial anisotropy modeled as Ising sys- 
tem which is from a numerical point of view an 
ideal candidate to study basic properties of EB. 
However, since the occurrence of EB is not re- 
stricted to systems with a strong anisotropy in 
the AFM, in the present paper we will extend the 
previous model @,|,|K} to the non-Ising case, i. e. 
we will vary the strength of the anisotropy of the 
AFM. 

In the next section we give a brief review of 
the physics of domains in diluted Ising antiferro- 
magnets in an external field. These systems have 
been studied in great detail in the past and the 
physics which emerge from these studies are im- 
portant for understanding EB. In Sec. || our model 
is described and in Sec. f| our results from Monte 
Carlo simulations are discussed. Finally, we sum- 
marize in the last section. 

2 Domains in disordered an- 
tiferromagnets 

Considerable interest has been focused in re- 
cent years on the understanding of diluted Ising 
antiferromagnets in an external magnetic field 
(DAFF) as they are ideal candidates for the study 
of disordered systems. A typical material for ex- 
perimental studies is Fei_ p Zn p F2 where FeF2 is 



the AFM which is randomly diluted with proba- 
bility p by non-magnetic Zn ions. Theoretically, 
due to the very strong uniaxial anisotropy this 
system is usually treated as Ising model. Prop- 
erties which have been extensively exploited are 
the critical behavior, domain structures, metasta- 
bility and slow dynamics (for reviews on DAFF 
see 00). Additionally, many of the findings of 
the DAFF are also relevant for the random field 
Ising model (RFIM) which has been shown to be 
in the same universality class 

In zero field the system undergoes a phase tran- 
sition from the paramagnetic phase to the long- 
range ordered, antiferromagnetic phase at the dis- 
order dependent Neel temperature Tjv as long as 
the dilution p is small enough so that the lattice 
of occupied sites is above the percolation thresh- 
old. In the low temperature region, for small mag- 
netic fields B the long-range ordered phase re- 
mains stable in three dimensions ||l5|, fig], while 
for higher fields the DAFF develops a disordered 
domain state 17, 18) with a spin-glass-like behav- 
ior. The reason for the domain formation was 
originally investigated by Imry and Ma for the 
RFIM p9| . The driving force for the domain for- 
mation is a statistical imbalance of the number of 
impurities of the two antiferromagnetic sublattices 
within a finite region of the DAFF. This leads to 
a net magnetization of this region which couples 
to the external field. A spin reversal of this re- 
gion, i. e. the creation of a domain can hence 
lower the energy of the system. The necessary en- 
ergy increase due to the formation of a domain 
wall can be minimized if the domain wall passes 
preferentially through non-magnetic defects at a 
minimum cost of exchange energy. Hence, these 
domains have non-trivial shapes following from 
an energy optimization. They have been shown 
to have a fractal structure with a broad distribu- 
tion of domain sizes and with scaling laws quan- 
titatively deviating from the original Imry-Ma as- 
sumptions [po| , ^i| . 

In small fields the equilibrium phase of the 
three-dimensional DAFF is long-range ordered. 
However, if cooled in a field B below a certain 
temperature Ti(B), the system usually also devel- 
ops metastable domains H). The reason for 
this metastability is a strong pinning which hin- 
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ders domain wall motion. These pinning effects 
are due to the dilution (random-bond pinning) as 
well as due to the fact that a rough domain wall 
also carries magnetization in a DAFF (following 
again the Imry-Ma argument) which couples to 
the external field and hinders domain wall mo- 
tion (random- field pinning) Ipifl , Consequently, 
after cooling the system from the paramagnetic 
phase within an external field, a DAFF freezes in a 
metastable domain state which survives even after 
switching off the field, then leading to a remanent 
magnetization which decays extremely slow 25 1. 

In the following we will argue that these well 
established properties of the DAFF are the key 
for understanding exchange bias. Indeed, during 
preparation of the system, the AFM is cooled in 
an external magnetic field and additionally un- 
der the influence of an effective interface exchange 
field stemming from the magnetized FM. Hence, 
the AFM — as far as it is diluted in any sense — 
must develop a domain state with a surplus mag- 
netization similar to that of a DAFF after field 
cooling. 



3 Model for exchange bias 

The model which we consider in the following con- 
sists of one FM monolayer exchange coupled to a 
diluted AFM film consisting of t monolayers. In 
Fig. [I] a sketch of our model is shown for t = 3. 




Figure 1: Sketch of the model with one FM layer 
and three diluted AFM layers. The dots mark 
defects (non-magnetic ions). 

The system is described by a classical Heisen- 
berg model, 



n 



■Jfm ~ S (dzSf z +d x Sf x +S i -B_j 



(hi) 



-Jafm e i e j°Li-Vj ~ (kz^iz + 

— JlNT £i<Li ■ S_j, 

(id) 

where the first line contains the energy contribu- 
tions of the FM. Here, the first term is the fer- 
romagnetic nearest neighbor interaction with ex- 
change constant Jfm- The second term intro- 
duces an easy axis in the FM (z-axis, anisotropy 
constant d z — O.IJfm) which sets the Stoner- 
Wohlfarth limit of the coercive field, i. e. the 
zero temperature limit for magnetization rever- 
sal by coherent rotation (B c = 2d z , in our units, 
for a field parallel to the easy axis). The shape 
anisotropy is approximated by the next term 
(anisotropy constant d x — —O.IJfm) leading to a 
magnetization which is preferentially in the y — z- 
plane. We checked, however, that its value does 
not influence our results. The last term of the first 
line is the Zeemann energy. 

The second line describes the diluted AFM 
(ei = 0,1; dilution p) correspondingly except of 
the shape anisotropy. For the exchange constant 
of the AFM which mainly determines its Nccl tem- 
perature (also depending on the dilution and the 
uniaxial anisotropy k z ) we set Jafm = — Jfm/2. 
Finally, the third line includes the interface cou- 
pling between FM and AFM and for simplicity we 
assume Jint = — Jafm)- 

Our magnetic field B will always be along the z 
axis. In earlier publications the AFM was 

described by an Ising model. In the present work, 
we relax this restriction on the AFM. In order to 
investigate a broader class of systems for the AFM 
we vary the uniaxial anisotropy k z of the AFM. 

4 Results from Simulations 

We use Monte Carlo methods with a heat-bath al- 
gorithm and single-spin flip methods for the sim- 
ulation of the model explained above. Each spin 
is subject to a trial step consisting of a small de- 
viation from the original direction followed by a 
second trial step in form of a total flip. This two- 
fold trial step can take care of a broad range of 
anisotropics starting from very soft spins up to 
the Ising limit |2(J. We perform typically 25000 
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Monte Carlo steps per spin for a complete hystere- 
sis loop. 

Since we are not interested in the critical be- 
havior of the model above, we do not perform a 
systematic finite-size analysis. However, in order 
to observe the domain structure of the AFM we 
have to guarantee that typical length scales of the 
domain structure fit into our system. Therefore, 
we show here only results for rather large systems 
of lateral extension 128 x 128. Nevertheless, we 
also varied the system size and checked that our 
results are not influenced by the system size as 
long as the system is not much smaller. 

In our simulations the system is cooled from 
above to below the ordering temperature of the 
AFM. During cooling the FM is long-range or- 
dered along the easy z axis and its magnetization 
is practically constant, resulting in a nearly con- 
stant exchange field for the AFM monolayer at 
the interface. In addition to this exchange field 
the external, magnetic field acts also on the whole 
AFM. 




Figure 2: Frozen domain states of a 40% di- 
luted AFM consisting of 6 monolayers for dif- 
ferent values of the AFM anisotropy, k z = 
0.1 Jfm, 1-0 Jfm, 30 Jfm (from top). The shading 
codes the z-component of the staggered magneti- 
zation. 



As already argued in the section before, during 
the cooling procedure the AFM becomes frozen in 
a domain state, the structure of which depends 
on the system parameters. The influence on di- 
lution |0| and the influence of the AFM film 
thickness was already discussed before for the 
case of an Ising AFM. In the present case, typ- 
ical staggered domain configurations of the bulk 
AFM are shown for three different values of the 
AFM anisotropy (Fig. ||). For low anisotropies 
k z < Jafm domain walls have a width of the or- 
der of yj J afm /kz ■ Even for the lowest anisotropy 
shown in Fig. [2| the width is only of the order of a 
few lattice constants which can hardly be detected 
in our figure. Also, due to the dilution walls tend 
to follow the holes so that the wall width is further 
reduced at those places. Interestingly, the domain 
structure itself depends also on k z . The system 
has the smallest domains for an intermediate value 
of k z = Jfm and not for the Ising case correspond- 
ing to the high anisotropy limit k z = 30 Jfm as one 
might expect. We will discuss the results follow- 
ing from this behavior later in connection with the 
anisotropy dependence of the EB. 

Typical hysteresis loops taken after cooling in a 
field of B = 0.25 Jfm are depicted in Fig. ||. Shown 
are results for the magnetization of the FM (up- 
per figure) as well as that of the AFM interface 
layer (lower figure). An exchange bias is observed 
clearly and we determine the corresponding ex- 
change bias field as Beb = {B + + B~)/2 where 
B + and B~ are those fields of the hysteresis loop 
branches for increasing and decreasing field, where 
the easy axis component of the magnetization of 
the FM becomes zero. 

The interface magnetization of the AFM also 
shows a hysteresis, following the coupling to the 
FM. Additionally, its curve is shifted upwards due 
to the fact that after field cooling the AFM is in a 
domain state with a surplus magnetization. The 
upward shift of the hysteresis loop for the interface 
AFM proves the existence of remanent magnetiza- 
tion in the AFM domains. This shifted interface 
magnetization of the AFM acts as an additional 
effective field on the FM, resulting in EB. The 
magnitude of the EB field strongly depends on 
the amount of this upward shift. Note that the 
shift is of the order of a few percent of the satu- 
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Figure 3: Typical hysteresis loop along z of a) the 
FM and b) the interface layer of the AFM, for a 
dilution p = 0.4 and an AFM thickness t = 2. The 
magnetization is in units of the saturation value 
and the field in units of Jfm- 



ration magnetization of the AFM while approxi- 
mately 10% of the spins of the AFM contribute to 
the AFM hysteresis. The saturation field for the 
AFM is much larger than that of the FM so that 
the AFM is never saturated during the simulation. 

In absence of any anisotropy in the FM and at 
very low dilution of the AFM we observe a per- 
pendicular coupling between FM and AFM. The 
magnetization reversal in the FM is here by coher- 
ent rotation. The picture changes with increasing 
uniaxial anisotropy in the FM and upon further di- 
lution of the AFM. The magnetization reversal in 
the FM is now by domain wall motion and the per- 
pendicular coupling becomes less significant. This 
is because uniaxial anisotropies in both the FM 
and the AFM having the same axis no longer lead 
to an energy minimum for a perpendicular cou- 
pling across the interface. Moreover AFM spins 
with missing AFM neighbors can lower their en- 
ergy by rotating parallel to their FM neighbors. 
Therefore, in the framework of our calculations a 
spin-flop coupling is not an essential mechanism 
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Figure 4: EB field versus anisotropy of the AFM 
for different AFM thicknesses (numbers of AFM 
layers). Field and anisotropy are in units of Jfm- 

We have calculated the EB field for a wide range 
of values of k z , starting from very soft spins to 
rigid, Ising-like spins. Fig. || shows result for dif- 
ferent thicknesses of the AFM and for a dilution 
of p = 0.4. Interestingly, we find a peak in the 
EB field at an intermediate value of k z for a suffi- 
ciently thick AFM while at lower thicknesses the 
EB field increases with the anisotropy and satu- 
rates in the Ising limit. 

The key for the understanding of EB is the un- 
derstanding of AFM domain configurations and 
domain walls. AFM domains are required to carry 
a surplus magnetization at the interface which 
must be stable along the z direction during hys- 
teresis in order to produce any EB. In general, one 
might expect that the most stable domain con- 
figurations are obtained for the Ising limit. But 
Fig. H and also Fig. || suggest that the behav- 
ior of domains is more complex. Let us start 
considering the Ising limit where some domain 
wall is formed upon field cooling. This domain 
wall preferentially passes through defects thereby 
minimizing the exchange energy and at the same 
time it gathers magnetization thereby lowering the 
Zeeman energy. When the anisotropy k z is de- 
creased the energy to create a domain wall will 
decrease. Thus the system will respond by rough- 
ening the domain boundaries (see Fig. |^, where 
the domain configuration in the middle is more 
complex than the lower one which represents the 
Ising limit). This roughening enhances the possi- 
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bility for the domains to carry any surplus magne- 
tization and hence the EB will increase. However, 
there exists a counter effect. While further de- 
creasing k z the width of the domain wall increases, 
so that less energy can be saved through the di- 
lution. Hence, for still lower anisotropy the do- 
main walls will smoothen thereby lowering again 
the exchange energy (see once again Fig. |2j, now 
comparing the domain configuration in the mid- 
dle with the upper one for still lower anisotropy 
which has much smoother domain walls). Since 
flat walls carry less remanent magnetization the 
EB will decrease now with decreasing anisotropy. 
The compromise between these two opposite ef- 
fects is achieved at some intermediate value of k z 
where the bias shows a peak. 

However, the peak disappears at lower values 
of t. This happens since for only one monolayer 
of AFM we are close to the percolation thresh- 
old where the domain walls pass nearly exclusively 
through the defects costing very little or no energy. 
Therefore the first mechanism discussed above is 
less important and the EB increases with k z till it 
saturates in the fsing limit. 



5 Conclusions 

In conclusion, we find that the domain state model 
for EB proposed originally for the Ising AFM is 
not restricted to this limit. Rather under cer- 
tain combination of thickness and dilution of the 
AFM, the softness of the AFM spins can lead to 
an even stronger bias field. Since disorder in the 
AFM of an exchange bias system is rather com- 
mon, our model yields a general understanding of 
the microscopic origin of exchange bias. Within 
our model there are several properties which influ- 
ence the bias field, such as dilution, thickness and 
anisotropy of the AFM. Although a qualitative un- 
derstanding regarding the dependence of EB on 
these parameters has been achieved, a quantita- 
tive study of the domain structure both at the in- 
terface and in the bulk of the AFM would provide 
a deeper understanding to the problem. 
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